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Abstract 

We investigate the set of spacetime general coordinate transformations (G.C.T.) 
which leave the line element of a generic Bianchi Type Geometry, quasi-form in- 
variant; i.e. preserve manifest spatial Homogeneity. We find that these G.C.T. 's, 
induce special time-dependent automorphic changes, on the spatial scale factor 
matrix *f a p(t) -along with corresponding changes on the lapse function N(t) and 
the shift vector N a (t). These changes, which are Bianchi Type dependent, form 
a group and are, in general, different from those induced by the group SAut(G) 
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advocated in earlier investigations as the relevant symmetry group; they are used 
to simplify the form of the line element -and thus simplify Einstein's equations 
as well-, without losing generality. As far as this simplification procedure is con- 
cerned, the transformations found, are proved to be essentialy unique. For the case 
of Bianchi Types II and V, where the most general solutions are known -Taub's 
and Joseph's, respectively-, it is explicitly verified that our transformations and 
only those, suffice to reduce the generic line element, to the previously known 
forms. It becomes thus possible, -for these Types- to give in closed form, the most 
general solution, containing all the necessary "gauge" freedom. 



1 Introduction 

It is well known, that spatial homogeneity reduces Einstein's field equations for pure 
gravity, to a system of ten coupled O.D.E.'s with respect to time |]: one equation 
quadratic in the velocities 7 aj g and algebraic in N 2 (G o = 0), three linear in velocities 
and also algebraic in N a (G 0i = 0), and the six spatial equations (Cry = 0) which are 
linear in j a/ 3 and are also involving N, N, N a , N a , j a p, ^ a p. 

In attempting to find solutions to this set of equations, it is natural -although seldom 
adopted in the literature- to solve the quadratic constraint for N 2 and the linear con- 
straint equations for as many of the iV a 's as possible; then substitute into the remaining 
spatial equations. When this is done, the spatial equations can be solved for only 6-4=2 
independent accelerations ^ a p. Only for Bianchi Type II and III -a particular VI case- 
we can solve for 6-3=3 accelerations, since only two of the three linear constraints are 
independent; but then in both of these linear combination of the N a, s remains 

arbitrary and counterbalances the existence of the third independent acceleration. Thus, 
the general solution to the above mentioned system of equations will, in every Bianchi 
Type, involve four arbitrary functions of time, whose specification should, somehow, cor- 
respond to a choice of time and space coordinates -in complete analogy to the full pure 
gravity, whereby four arbitrary functions of the spacetime coordinates, are expected to 
enter the general solution. 

In the literature a different approach is more frequently met. It involves an a priori gauge 
choice of coordinate system: As far as time is concerned, one may set N to be either an 
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explicit function of time -say 1 or t 2 e.t.c.-, or some combination of 7 a/ g's -see (2.8). For 
the spatial coordinates, the depicted situation is more vague. In some works, iV a, s are 
set to zero, in others, some N a 's are retained. In any case, most of these choices, are 
considered as being, more or less, inequivalent and their connection to the well-known 
existence of Gauss-normal coordinates (goo — — 1> 9oi — 0) 0! is n °t at all clear. When 
such a gauge choice has been made, the spatial equations can be solved for all 6 inde- 
pendent j a p{t)- The constraint equations become then algebraic equations, restricting 
the initial data -needed to specify a particular solution of the spatial equations. 
In both approaches, the ensuing equations are still too difficult to handle; thus further 
simplifying hypotheses are employed, such as N a (t) = 0, leading to 7 a/ 3 = diag((a 2 (t), 
b 2 (t) , c 2 (t))) for Class A Types e.t.c. For the Bianchi types I and IX, the hypothesis 
N a (t) = and j a p = diag((a 2 (t),b 2 (t),c 2 (t))), is known to be linked to kinematics 
and/or dynamics -although in a, somewhat, vague way see e.g. @ and Ryan in 
In all other cases, this or any other simplifying hypothesis used, is interpreted only as 
an ansatz to be tested at the end, i.e. after having solved all the -further simplified- 
equations. For example, to take an extreme case, diagonality of 7 a /j(t) together with the 
vanishing of the shift vector is known to lead to incompatibility for Bianchi types IV, 
VII (Class B) [g, |[, as well as for the biaxial type VIII cases (a 2 , a 2 , c 2 ), (a 2 , b 2 , a 2 ), |J. 
The diversity of the various ansatzen appearing in the literature, causes a considerable 
degree of fragmentation. 

It has long been suspected and/or known, that automorphisms, ought to play an 
important role in a unified treatment of this problem. The first mention, goes back 
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to the first of 0. More recently, Jantzen, -second of ||- has used Time Dependent 
Automorphism Matrices, as a convenient parametrization of a general positive definite 
3x3 scale factor matrix 7 Q!( a(i), in terms of a -desired- diagonal matrix. His approach, is 
based on the orthonormal frame bundle formalism, and the main conclusion is (third of 
[6], pp. 1138): "... the special automorphism matrix group SAut(G), is the symmetry 
group of the ordinary differential equations, satisfied by the metric matrix 7 a/ g ; when no 
sources are present ..." Later on, Samuel and Ashtekar in |7j , have seen automorphisms, 
as a result of general coordinate transformations. Their spacetime point of view, has led 
them, to consider the -so called- " Homogeneity Preserving Diffeomorphisms" , and link 
them, to topological considerations. 

In this paper, we also take a spacetime point of view, and try to avoid the frag- 
mentation -mentioned above-, by revealing those G.C.T.'s, which enable us to simplify 
the line-element -and therefore Einstein's equations-, while at the same time, they pre- 
serve manifest spatial homogeneity. We are, thus, able to identify special automorphic 
transformations of 7 a/ 3(t), along with corresponding -non tensorial, for the shift vector- 
changes of N,N a which allow us to set N a = and bring 7 a/ g(i) to some irreducible, 
simple -though not unique- form. 

The structure of the paper, is organized as follows: 
In section 2, after establishing the existence of Time-Dependent Automorphism Induc- 
ing Diffeomorphisms (A.I.D.'s), the general irreducible form of the line element for all 
Bianchi Types is given, and a uniqueness theorem, is proven. 

In section 3, attention is focused on Bianchi Types II and V. Einstein's equations obtain- 
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ing from the irreducible form of the line element, are explicitly written down and com- 
pletely integrated. The uniqueness of the transformations given in section 2, is explicitly 
verified, with the aid of the well known Taub's and Joseph's solution -respectively. As 
a result, we give the closed form of the most general line elements, satisfying equations 
(2.5). 

Finally, some concluding remarks are included in the discussion. 

2 Time Dependent Automorphism 
Inducing Diffeomorphisms 

It is well known that the vacuum Einstein field equations can be derived from an action 
principle: 

-1 



A = Ttr V-^i^Rd'x (2.i) 

107T J 

(we use geometrized units i.e. G = c = 1) 

The standard canonical formalism || makes use of the lapse and shift functions appearing 
in the 4- metric: 

ds 2 = (N*Ni - N 2 )dt 2 + 2N i dx i dt + g ij dx i dx j (2.2) 

From this line-element the following set of equations obtains, expressed in terms of the 
extrinsic curvature: 



Ho = Jg~ {KijKV -K 2 + R) = 



(2.3a) 



Hi = (K^ - K\i) = (2.3b) 

-^=f t iV9 - Kg*)] = -N{R» - l -R g*) - j(K kl K kl - K 2 )g* 

+2N(K ik K{ - K K ij ) - (N lij - N^g*') + [(K ij - K g ij )N% ( 2 - 3c ) 

-NfiiK 1 * - K g lj ) - N^(K li - K g li ) 
This set is equivalent to the ten Einstein's equations. 

In cosmology, we are interested in the class of spatially homogeneous spacetimes, 
characterized by the existence of an m-dimensional isometry group of motions G, acting 
transitively on each surface of simultaneity E t . When m is greater than 3 and there is 
no proper invariant subgroup of dimension 3, the spacetime is of the Kantowski-Sachs 
type and will not concern us further. When m equals the dimension of E 4 -which is 
3-, there exist 3 basis one-forms cr" satisfying: 

do" = C| 7 /A.^ a? j -o? ti = 2C% aj a? (2.4a) 

where are the structure constants of the corresponding isometry group. 

In this case there are local coordinates t, x l such that the line element in (2.2) assumes 

the form: 

ds 2 = (N a (t)N a (t) - N 2 (t))dt 2 + 2N a (t)o?(x)dx i dt 

(2.4b) 

+ r ) a f3(?i{x)vj{x)dx l dx j 
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Latin indices, are spatial with range from 1 to 3. Greek indices, number the different 
basis 1-forms, take values in the same range, and are lowered and raised by j a p, and 7 a/3 
respectively. 

A commitment concerting the topology of the 3-surface, is pertinent here, especially 
in view of the fact that we wish to consider diffeomorphisms 0; we thus assume that 
G is simply connected and the 3-surface Sf can be identified with G, by singling out a 
point p of S t , as the identity e, of G. 

If we insert relations (2.4) into equations (2.3), we get the following set of ordinary 
differential equations for the Bianchi-Type spatially homogeneous spacetimes: 

E = K% Kl - K 2 + R = (2.5a) 



E a = K£ C; e - K* C e afl = (2.5b) 

E% = K% - NKKp + NR« + 2N"{K - K» C a up ) (2.5c) 
where K* = r y ap K p p and 

K aP = —{iap + 2 lau C^ p NP + 2 lpv Cl p N") (2.6) 



R a p = C^C^y^l^ + ZCifipx + ZC&CfrYrT* 

When N a = 0, equation (2.5c) reduces to the form of the equation given in jlOfl . Equation 
set (2.5), forms what is known as a -complete- perfect ideal; that is, there are no 



(2.7) 



integrability conditions obtained from this system. So, with the help of (2.5c), (2.6), 
(2.7), it can explicitly be shown, that the time derivatives of (2.5a) and (2.5b) vanish 
identically. The calculation is staightforward -although somewhat lengthy. It makes use 
of the Jacobi identity CfyC^+C^C^+C^C^ = 0, and its contracted form C^C^ = 0. 

The vanishing of the derivatives of the 4 constrained equations: 
E = 0,E a = 0, implies that these equations, are first integrals of equations (2.5c) 
-moreover, with vanishing integration constants. Indeed, algebraically solving (2.5a), 
(2.5b) for N(t),N a (t), respectively and substituting in (2.5c), one finds that in all -but 
Type II and III- Bianchi Types, equations (2.5c), can be solved for only 2 of the 6 
accelerations 7 Q/ g present. In Type II and III, the independent accelerations are 3, since 
E a are not independent and thus, can be solved for only 2 of the 3 N a, s. But then in both 
of these linear combination of the iV a 's remains arbitrary, and counterbalances 

the extra independent acceleration. Thus, in all Bianchi Types, 4 arbitrary functions 
of time enter the general solution to the set of equations (2.5). Based on the intuition 
gained from the full theory, one could expect this fact to be a reflection of the only known 
covariance of the theory; i.e. of the freedom to make arbitrary changes of the time and 
space coordinates. 

The rest of this section is devoted to the investigation of the existence, uniqueness, 
and properties of general coordinate transformations -containing 4 arbitrary functions 
of time-, which on the one hand, must preserve the manifest spatial homogeneity, of the 
line element (2.4b), and on the other hand, must be symmetries of equations (2.5). 
As far as time reparametrization is concerned the situation is pretty clear: If a transfor- 
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mation 



t^t = g(t)&t = /(f) 



(2.8a) 



is inserted in the line element (2.4b), it is easily inferred that 



(2.8b) 




dt 



(2.8c) 



Accordingly, Kp transforms under (2.8a) as a scalar and thus (2.5a), (2.5b) are also 
scalar equations while (2.5c) gets multiplied by a factor df(t)/dt. Thus, given a par- 
ticular solution to equations (2.5), one can always obtain an equivalent solution, by 
arbitrarily redefining time. Hence, we understand the existence of one arbitrary function 
of time in the general solution to Einstein's equations (2.5). In order to understand 
the presence of the rest 3 arbitrary functions of time it is natural to turn our atten- 
tion to the tranformations of the 3 spatial coordinates x\ To begin with, consider the 
transformation: 



It is here understood, that our previous assumption concerning the topology of G and 
the identification of St with G, is valid for all values of the parameter t, for which the 
transformation is to be well defined. 



t = t&t = t 



(2.9) 



x l = g\x j ,t) <^x* = f{x j ,t) 



8 



(2.10) 



Under these transformations, the line element (2.4b) becomes: 

ds 2 = [(N a N a - N 2 ) + of (/)7«*(t) 

+ 2a?U)%N a (t)}dt 2 

Since our aim, is to retain manifest spatial homogeneity of the line element (2.4b), we 
have to refer the form of the line element given in (2.10) to the old basis of (x) at the 
new spatial point x % . Since erf -both at x l and well as, dx l /dx\ are invert ible 

matrices, there always exists a non-singular matrix A*(x,t) and a triplet P a (x,t), such 
that: 

°^ . (2.11) 

f(x)^=P a (x,t) 
ot 

The above relations, must be regarded as definitions, for the matrix A° and the triplet 
P a . With these identifications the line element (2.10) assumes the form: 

ds 2 = [(N a N a - N 2 ) + P a (x,t)pP(x,t)^ a p(t) + 2P a (x,t)N a (t)]dt 2 

+ 2A°(x, t)a^(5;)[N a (t) + pP(x, t)^(t)}dx m dt (2-12) 
+ A«(5, t)A^(x, t)^(t)aUx)<(x)dx m dx n 
If, following the spirit of [0, we wish the transformation (2.9) to be manifest homogeneity 
preserving i.e. to have a well defined, non-trivial action on j a p{t), Af(t) and N a (t), we 
must impose the condition that A™ (a;, t) and P a (x, t) do not depend on the spatial point 
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x, i.e. A° = A°(t) and P a = P a (t). Then (2.12) is written as: 

ds 2 = [(N a N a - N 2 ) + P a P^ a/3 + 2P a N a }dt 2 
+ 2A^(x)[N a + pP laP ]dx m dt 

+ A a Xlafs<(x)a:(x)dx m dx n (2-13) 
ds 2 = (N a N a - N 2 )dt 2 + 2N a (t)a?(x)dx i dt 

+ lap{t)cr^ (x)(Tj (x)dx l dx j 

with the allocations: 

7^ = A^7m, (2.14a) 
N a = A^Np + P P l P n) and thus N a = S^(N^ + P p ) (2.14b) 

N = N (2.14c) 

(where S = A' 1 ). 

Of course, the demand that and P a must not depend on the spatial point x\ changes 
the character of (2.11), from identities, to the following set of differential restrictions on 
the functions defining the transformation: 

^ = <U)m°^) (2-15a) 

% = <U)P a {t) (2.15b) 

Equations (2.15) constitute a set of first-order highly non-linear P.D.E.'s for the un- 
known functions f\ The existence of local solutions to these equations is guaranteed by 
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Frobenius theorem Oj as long as the necessary and sufficient conditions: 



dxi\dx m ) dx m \dxi) 



d f dp \ d /dp 
dt\dx m ) dx m \dt 











hold. Through (2.15) and repeated use of (2.4a), these equations reduce respectively to: 

K c h = K K C % ( 2 - 16 ) 



P*C°J% = \^ (2-17) 



It is noteworthy that the solutions to (2.16) and (2.17), -by virtue of (2.14)- form a 
group, with composition law: 

(Aa)2 = (Ai)^(A 2 )S 

(^) a = (Ai) a ,( p ^ + W 

where (A 1 , Pi) and (A 2 , P2), are two successive transformations of the form (2.14). 
Note also, that a constant automorphism is always a solution of (2.16), (2.17); indeed, 
Ajg(£) = A^ and P a (t) = solve these equations. Thus, A^ and P a = can be regarded 
as the remaining gauge symmetry, after one has fully used the arbitrary functions of time, 
appearing in a solution A%(t) and P a (t). Consequently one can, at first sight, regard all 
the arbitrary constants encountered when integrating (2.17), as absorbable in the shift, 
since the transformation law for the shift, is then tensorial. This is certainly true, as long 
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as there is a non zero initial shift. However, if one has used the independent functions 
of time, in order to set the shift zero, then the constants remaining within Ajg, are not 
absorbable. It is this kind of constants that we explicitly present below, when we give the 
solutions to (2.16), (2.17) for all Bianchi Types. A relevant nice discussion, distinguishing 
between genuine gauge symmetries (cf. arbitrary functions of time) and rigid symmetries 



cf. arbitrary constants), is presented in JT^|. There a different definition of manifest 



homogeneity preserving diffeomorphisms -stronger than the one adopted in this work- is 
used, and results in only the inner automorphisms being allowed to acquire t dependence. 
In connection to this, it is interesting to observe that (2.16-17) give essentially the same 
results: notice that 2P fl C® j3 is, by definition, the generator of Inner Automorphisms. 
Thus there is always a \%(t) = Exp{2P»C* p ) G IAut(G) satisfying (2.17). If we now 
parameterize the general solution to (2.16-17) by Ap(t) = A°(t)Z7|(t) and substitute in 
these relations, we deduce that the matrix U is a constant automorphism. This analysis 
is verified in the explicit solutions to (2.16-17), presented latter. 

Equation (2.16) is satisfied if and only if, A2(i) is an element of the automorphism 
group of the Lie algebra determined by the . Equation (2.17) further restricts the 
form of Aa(t) and P a (t), so that manifest spatial homogeneity is preserved despite the 
mixing of the old time and space coordinates in the new spatial coordinates x l . Thus, it 
is appropriate to call transformations (2.9), satisfying conditions (2.15), (2.16), (2.17), 
Time-Dependent Automorphism Inducing Diffeomorphisms. The importance of 
automorphisms in Bianchi Cosmologies, has been stressed in ||. The symmetry group of 
the differential equations, satisfied by 7^ (t), -advocated in these works of Jantzen et al- 
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is the unimodular matrices SAut(G). As we shall later see, we find another symmetry 
group, whose time-dependent part lies essentially in IAut(G) and thus coincides with 
SAut(G), only for Class A Bianchi Types VI, VII, VIII, IX. 

At this point it is natural to ask how this difference occurs. It is our opinion that the 
key elements on which the difference in the symmetry groups found rests, are: 

a) The inhomogeneous transformation law (2.14b) for the shift. Indeed, Jantzen 
(1979), having adopted an orthonormal frame-bundle point of view, naturally as- 
sumes as his "gauge" transformation laws (2.14a,c) and the tensorial law 

N a = S a p N p , for the shift p. 221. 

b) The different definition and/or role reserved for the triplet P a {t); we define it as 
a sort of "velocity" of the transformation (2.9) in (2.15b) and use it in the inho- 
mogeneous law (2.14b). On the other hand, Jantzen (1979), (p. 221) uses the 
corresponding quantity uj a {t) (so called velocity of the automorphism frame)to de- 
fine a new time derivative d/di = d/dt + uj a {t)a l a (x)d/dx l . 

c) We concentrate on the symmetries of the O.D.E.s (2.5), i.e. of Eistein's equations 
written in the invariant base, while Jantzen, as far as we understand, focuses on the 
symmetries of the P.D.E.s (2.3), i.e. of Einstein's equations, written in an arbitrary 
frame. 

In [0, the so called Homogeneity Preserving Diffeomorphisms, are considered in re- 
lation to the topology of S t . A time-independent version of (2.15), appears in [ |T^] . 



where all homogeneous three-geometries, are classified in equivalence classes, with re- 
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spect to these "frozen" transformations. It is straightforward to check, that E , E a , E% 
transform -under (2.14)- as follows: 

E = E , E a = A^Ep, E« = S^E2 (2.18) 

This fact, establishes the covariance of equations (2.5), under the "gauge" transformation 
(2.14), and implies that if (N, N a , 7^) is a solution to Einstein's equations, so will be 
the set (N, N a , j a/3 ) -provided that, (2.16), (2.17) hold-; in fact, as the preceding 
exposition proves, they will be the same equations expressed in different space-time 
coordinate systems. Out of the twelve quantities A^j(t) and P a (t), conditions (2.16), 
(2.17) leave us, as we are going to see, in every Bianchi Type, with 3 arbitrary functions 
of time. This fact, along with the time reparametrization covariance, completes our 
understanding of why four arbitrary functions of time enter the general solution to (2.5). 
Consequently, transformation (2.14), gives us the possibility to simplify the form of the 
line element, and thus of Einstein's equations without loss of generality. It is obvious, 
that the simplification obtained, is different for different Bianchi Types, and even within 
a particular Bianchi Type it is not unique -since one may "spend" the freedom of the 
three arbitrary functions in different ways. 

A particularly interesting result, is that the shift vector N a can always be put to 
zero -perhaps at the expense of a more complicated "f a p- For the sake of completeness, 
we give below, a detailed analysis of the space of solutions to (2.16) and (2.17), for each 



and every Bianchi Type (solutions to (2.16), have been presented in [14 



To this end, recall that in 3 dimensions, the tensor , admits a unique decomposition 
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in terms of a contravariant symmetric tensor density of weight -1, m a/3 and a covariant 
vector u a = \C^ p as follows |fi5| : 



C% = m a5 e 6M + Vfi 6° - vj% 

The contracted Jacobi identities imply that m a "vp = 0, i.e. u a is a null eigenvector of 
the matrix m al3 . Under a GL(3, 3?) linear mixing of the basis 1-forms a a — > a a = S^a 13 , 
the structure constant tensor transforms as: 



Accordingly, the m al3 and i/ Q transform as: 

A (and thus S 1 ) is called a Lie algebra automorphism if Cg = C^ 7 , i.e. if m a/3 and Z/ a are 
equal to m al3 and v a respectively. In this case the automorphism conditions become: 



m 



P = \S\~ 1 S^S^m' /S (2.19a) 



v a = A^p (2.19b) 

The different Bianchi Types, arise according to the rank and signature of m al3 , and 
the vanishing or not, of v a . Using (2.19), one can -straightforwardly- solve the system 
of equations (2.16) and (2.17). We now present, the form of A2(i) and P a (t) satisfying 
(2.16), (2.17) as well as some irreducible form for j a p, for each Bianchi Type: 
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Type I: m al3 — 0, v a — 0. This Type has been exhaustively treated, in the literature 
(HI) 0)- We only note that -since all are zero- (2.17), implies that P a (t) is arbitrary 
and Ajg(t) is constant. Then, (2.16) implies that is an element of GL(3, 9ft). Thus, 
without loss of generality, one can set iV Q = 0, -using (2.14b). A first integral of 
equations (2.5c) is then, 7 ap 7 Pj g = where is an arbitrary constant matrix. From 
this point, the standard textbooks, || deduce (using algebraic arguments) a diagonal 
form: 

la/3 — diag(e at , e' 3 *, e 7 *) and then using Einstein's equations find the general solution, 
which depends on 1 essential parameter, as expected -see table. 

Indeed, from (2.5c), one has 12 initial constants; 6 7 aj g, and 6 7 Q/ 3 at some t -according 
to Peano's theorem. The quadratic constraint equation (2.5a), reduces them to 10, and 
then, with the usage of constant automorphisms -which contain 9 A2's-, the number of 
the remaining essential constants (or essential parameters), is 10 — 9 = 1. 

Type II rank(m) = l and v a = . Then, matrix m al3 , can be cast to the form 



ni 



a/3 



diag(l/2,0,0). Equations (2.16), (2.17) give the following form for A*g(t): 



QiQa ~ Q203 x(t) y(t) 



AS(«) 







^1 Q2 



\ 



Q3 Q4: 



[Qi, 62, 63, Q4 constants) 



The triplet P a (t) assumes the form: 



6164-6263 6164-6263 



The general solution to this Type, is Taub's solution (p6|), which contains 2 essential 
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parameters -see table. 



Again, we can understand this number, using Peano's theorem and the arbitrary extra 
constants, appearing in Ajg. Using x(t) and y(t), we start with 4 7^'s (i.e. we set 
712 = 7i3 = 0) and no shift. Thus the initial arbitrary constants, are 2x4 = 8. Out 
of these, the quadratic constraint equation (2.5a), removes 2, and 4 more are eliminated 
by the 4 g : s, contained in Ap 1 . So, the remaining arbitrary constants are: 8 — 2 — 4 = 2, 
in accordance with the number of expected essential parameters. 

Type V rank(m)=0 and v a ^ 0. Then m al3 = and we may choose v a = — hS% 
Equations (2.16), (2.17) give the following form for Ap(t): 



2 a' 



1 QiP(t) g 2 P{t) x(t) N; 



Q 3 P(t) QiP(t) y(t) 
1 J 

with Q1Q4 — Q2Q3 = 1 and the triplet: 



V 



'f?i) Q2, Q3, Qi constants) 



P a (t) = (x(ln|y, y(ln|)-, (lnl)) 



The general solution, is also known, as Joseph's solution (113)) with one essential pa- 
rameter. 

This number comes naturally, within our method; using x(t) and y{t), one can eliminate 
713 and 7 23 . Then, P(t) can serve to set the sub determinant of j a p, equal to (73a) 2 . At 
this stage, we are left with 3 r y a /3 ,s while the linear constraints equations (2.5b), imply 
that the shift is zero. Again, the quadratic constraint (2.5a), subtracts 2 arbitrary con- 
stants, and the constants contained in A^, 3 more. Then, the result is: 6 — 2 — 3 = 1, 
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essential constant. 



Type IV rank(m)=l and v a 7^ 0. We may choose 



m' 



a/3 



my 



diag(l/2,0,0), v a = -\5l- Equations (2.16), (2.17) give the following form for 



^ P{t) P{t) ln[/cP(t)] x(t) * 
k a {t) = o P(t) y(t) , ( K constant) 

1 J 



and the triplet: 



P"(t) = (x(\n^y-y, y(l n |)', (ml)') 



In this Type -which is a class B Type-, we can set 713 = 723 = 0, using x(t) and y(t). 
At this stage, the 2 of the 3 linear constraint equations, imply iV 1 = iV 2 = 0, while 
the third, involves P(t). Thus we can further, either set iV 3 = -through (2.14b)- and 
retain a non-zero 712, or eliminate 712, at the expense of a non- vanishing iV 3 . It is well 
known, that iV 3 = and 712 = 0, leads to incompatibility [4]. 
We have thus, the following counting of the essential parameters: 

a) When 7 12 7^ and N 3 = 0, we have 8 — 2 (from the quadratic constraint) — 2 (from 
the remaining linear equation) — 1 (from the constant contained in A^) = 3. 

b) When 712 = and iV 3 7^ 0, we have 6 — 2 (from the quadratic constraint) — 1 (from 
the constant contained in A^j) = 3. Notice that here, the remaining linear constraint 
equation, simply serves to define A^ 3 and thus, does not remove any constant. 



Type VI (Including Type III ||19|| , fT8|l ) rank(m)=2, signature (m)=Lorentzian and 
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v a 7^ 0. One convenient choice is m Q/3 = diag(l, —1,0) and v a = hS^. 
Note: In the standard texts e.g. /|J5|/, the matrix m a/3 is given in a more complicated 
form, which carries part of the arbitrariness of the magnitude of the vector v a . In this 
work, we imply that are given by their defining relation in terms of £ a/ 3 7 , m a/3 , v a . 
For all values of h ^ 0, ±1, equations (2.16), (2.17) give the following form for A^(t): 



AS(f) 



^ e~ hP ^Xcosh(P(t)) e- hp WAsinh(P(t)) x(t) X 
e- hP ^Xsmh{P{t)) e- hP ^\cosh{P(t)) y{t) 
1 



(A constant) 



while the triplet: 



P a (t) =(- 



(fr 2 - l)x(t)P(t) + fex(t) + y(t) 
2{h 2 - 1) ' 
(/i 2 -l)j/(*)P(0 + /iy(t)+i(t) P(t). 



2(fr 2 - 1) 

For h = 0, -class A-, there are two solutions: 

/ 



Acosh(P(t)) Asinh(P(t)) x{t) 
eAsinh(P(t)) eAcosh(P(t)) y(t) 



\ 



V 



o 







(A constant) 



while the triplet: 



P a {t) 



,ey(t) - x(t)P(t) ex(t) - y(t)P{t) eP(t) 
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where e = ±1. 



For h — ±1, -class B , the solutions are: 



f e -w(t)Acosh(P(i)) e- AP WAsinh(P(t)) 

e -^(*)Asinh(P(t)) e- ftp WAcosh(P(t)) c - hx(t) 



(A constant) 

while the triplet: 

_ 2MXg - cP(t) + 2hx(t)P(t) + ^ 

e-^^A sinh(P(t))P(f) - he~ hP ^X cosh(P(t))P(t) 
2e-^WAcosh(P(t)) - 2/ i e- ftp WAsinh(P(t)) ^ 

For each of the previously mentioned cases, we have: 

a) When h = 0, (class A), 7^ can be made diagonal and then the shift vanishes. Thus 
the counting of the essential parameters is: 6 — 2 (from the quadratic constraint) — 1 
(from the constant, contained in Ajg) = 3. 

b) When h = ±1, (class B), using x(t) and P(£), we can eliminate 713 and 7 23 . So: 8 
— 2 (from the quadratic constraint) — 2 (from the constants, contained in A^j) = 4 is 
the number of the essential constants. Notice that the 3 linear constraint equations, are 
linearly dependent and thus, when iV 3 = through (2.14b), there is no linear constraint 
equation left, to remove any constants, hence the number 4. 

c) When h ^ 0, ±1, the counting algorithm is exactly the same, as in Type IV case. 

Type VII rank(m)=2, signature (m)=Euclidean and u a 7^ 0. We may set m a/3 = 
diag{—\, —1, 0), v a — M 3 . For all values of h, equations (2.16), (2.17) give the following 

20 



form for A2(i): 



^ Xe hP ^ cos(P(t)) \e hP ® sin(P(t)) x(t) ^ 
_A e w(*) sin(P(0) Ae ftp W cos(P(t)) j/(t) 



1 y 



(A constant) 



and the triplet: 



P a {t) =( 



x(t)P(t) + h 2 x(t)P(t) - hx(t) + 



2(1 + /i 2 ) 

y(t)P(t) + h 2 y(t)P(t) - hy(t) ~ i(t) Pit) 
2(1 + /i 2 ) ' 2 J 

For the case h — 0, there is another solution, except the one deduced from the 



previous, by setting h — 0: 



Ag(*) 



Acos(P(t)) Asin(P(t)) 
Asin(P(t)) -Acos(P(t)) 

0-1 



(A constant) 



and the triplet: 



P a {t) 



,x(t)P(t)-y(t) y(t)P(t)+x(t) P(t) 



) 



Again, for each of the previously mentioned cases, we have: 
a) When h — 0, (class A), j a/3 can be made diagonal and equations (2.5b) give N a = 0. 
Thus: 6 — 2 (from the quadratic constraint) — 1 (from the constant, contained in A^) 
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= 3 is the number of the essential constants. 

b) When h ^ 0, the counting algorithm is exactly the same, as in Type IV case. 

For Bianchi Types VIII and IX, condition (2.17), does not impose any restriction on 
Ap(t), but rather fixes completely, the triplet P a {t) , to be: 

P a = -^—e^ TK m a ' 3 A T K^ s 
4|m| ' 

Type VIII rank(m)=3, signature (m)=Lorentzian. A standard choice is m Q/3 = 
i] al3 = diag(—l, 1, 1). Since \m al3 \ = — 1, (2.19a) implies that |A^| = 1 and thus, A^'s are 
the isometries of the Minkowski metric, in three dimensions, i.e. the Lorentz boosts, with 
one timelike and two spacelike directions, times a rotation of the "space" plane. Thus, 
the automorphisms are characterized by the two components of the velocity vector, plus 
the rotation angle. The triplet P a , is: 

P a = ^(AjA^r , -AJAJiT, -AjA^iT) 

It can be proven -see appendix A- that a positive definite matrix, can be diagonalized 
by this automorphism group; i.e. we can set 

7 Q/3 = diag(a 2 (t),b 2 (t),c 2 (t)). Then, from (2.5b), we will have N a = 0. 

The number 4, of the expected essential parameters -see table below-, can be under- 
stood as follows: The time-dependent Lorentz transformation A|g, can diagonalize j a ^. 
Thus, the counting: 6 — 2 (from the quadratic constraint) = 4. 

Type IX rank(m)=3, signature (m)=Euclidean. The standard choice is m alS = 5 al3 . 
Since |m a/3 | = 1,(2. 19a) implies that |A^| = 1 and thus, A^'s are the isometries of the Eu- 
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clidean metric, in three dimensions, i.e. the orthogonal matrices, which are characterized 
by three parameters; e.g. the Euler angles. The triplet P a , is: 

Since a positive definite symmetric matrix can be diagonalized by an element of -the con- 
nected to the identity component of- 0(3), we have that ■japit) = diag(a 2 (t), b 2 (t), c 2 (t)) 
19| . Then, from (2.5b), as is well-known N a = 0. 

The counting yields -exactly as in Type VIII: 6 — 2 (from the quadratic constraint) 
= 4, essential constants. 

^From the above analysis of the space of solutions to (2.16) and (2.17), we observe 
that in each Bianchi Type, there are 3 arbitrary functions of time -as expected-, for a 
twofold reason; 

Firstly, because we are solving the integrability conditions for the existence of a time- 
dependent spatial diffeomorphism according to (2.9). 

Secondly, because as it has been mentioned in the Introduction, the system of Einstein's 
equations (2.5), has a gauge freedom of 4 arbitrary functions of time. But one of them, 
simply corresponds to time reparametrization, while the remaining 3, are the ones we 
found in the above analysis. 

In the various Bianchi Types, the 3 arbitrary functions, are distributed differently among 
the components of A^(t) and P a (t). This fact, together with the different number of 
arbitrary constants appearing in A^ for each Type, results in a different number of 



essential constants -expected by independent arguments [19| to appear in the general 
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solutions to Einstein's equations (2.5) -see Table. 

We now conclude section 2, by stating the following (uniqueness) Theorem: 
"In a given -albeit arbitrary- Bianchi Type, let 7 1; j 2 , (in matrix notation) be solutions 
to Einstein's equations (2.5), then there is a matrix M of the form: M = A]" 1 EA 2 (where 
A x and A 2 are solutions to (2.16) and (2.17) and S 7 represents the irrelevant symmetries 
of the solution in its irreducible form) which connects them as: 72 = \ f ~ \\l ." 

Note: N, N a , are understood to be given from the quadratic and linear constraint 
equations (2.5a,b). 

The proof rests on the previously established facts: 

a) That the solutions to (2.16) and (2.17), suffice to reduce the generic 7 a/ j, to a form that 
will contain the expected necessary number of essential constants, so as to be regarded 
as the most general one -for each and every Bianchi Type. 

b) That for every given Bianchi Type, the solutions to (2.16) and (2.17), form a group. 
Indeed, let 71, 72 be solutions to (2.5). Then there are Ai, A 2 -along with Pi, P 2 
respectively, if needed- solutions to (2.16) and (2.17), such that: 

Tl '") 1 irreducible^-! 
72 = A^'Jirreducible^ 

where ^irreducible, stands for the solution in a form exhibiting, only the essential constants. 
From the first of these: 

^irreducible (Aj ) 71^-1 
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Since, -by definition- ^reducible is a symmetric matrix there are always, non-trivial 
matrices S, such that: 

1 irreducible ^ 1 irreducible^- 1 

Substituting the two last in the second, we obtain: 

72 = (A^ 1 SA 2 ) T 7l A 1 ' 1 SA 2 

q.e.d. 

3 The Space of Solutions for Type II and V Cases 

In this section, we adopt the more conventional point of view; that of "gauge" fixing, 
before solving. As far as time is concerned, we adopt the "gauge" fixing condition 
N = since this simplifies the form of the equations. For the spatial coordinates, as 
explained in the previous section, a choice of reference system, amounts to a choice of 
time-dependent automorphism -along with a choice of P a (t)-; thus, we select the form 
of l a p{t) to be such that, the linear equation would imply N a = 0. In this "gauge", 
Einstein's equations (2.5) read: 

_~-~/3A~ kA ~^ + ( 7 )2 _ m = Q (31a) 

Q,7 W 7 P£ - C; e Y P J pa = (3.1b) 
l af s - 7^7^7/9, - 27^/3 = (3.1c) 
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Note that taking the trace of equations (3.1c), one arrives at: 

(Zy-2yR=0 (3.2) 

7 

A somewhat useful result deriving from (3.2), is the following: 7 = ae 131 implies R = 0, 
which is incompatible for all but I Bianchi Types. 

We now present, a realization, of the method developed in the previous section, for the 
cases of Type II and V, Bianchi geometries. At this point, a word of caution is pertinent: 
it is evident -from the previously mentioned counting, of the expected number of essential 
constants-, that the well known Taub's (Type II) and Joseph's (Type V) solutions, are 
the most general for the respective cases [[LJ|. Thus, we should not expect to find 
anything new -in this respect. However, the thorough investigation of the complete 
space of solutions, requires the knowledge of the correct (gauge) symmetry group for 
Einstein's equations (2.5). In this respect, we shall directly show, that transformations 
(2.14), -as specified by conditions (2.16) and (2.16), applied to Types II and V-, are 
essentialy, the only (gauge) symmetries of these Bianchi geometries. 

Note: From now on we drop the tildes from the various quantities for simplicity 
-except in some cases, where misunderstanding could occur. 
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3.1 Bianchi Type II 



As it can be seen, from the results concerning Type II, we can consider -without loss of 
generality-, the time-dependent part 7^/3 (£), of the 3- metric, to have the form: 



a(t) 



b(t) f(t) 
\ /(f) c(t) J 

It is interesting to observe that, the freedom in arbitrary functions of time -contained in 
Ap(t)-, does not suffice to diagonalize 7 Q( g(£), i.e. to set fit) = 0, a priori. Yet, we know 
-see (3.16) and (3.17) below- that the diagonal Taub's metric, is the irreducible form of 
the most general Type II, solution. The reconciliation of these two, seemingly conflicting 
facts, obtains only on mass shell; after we have completely solved (3.1), with 7 Q( g(£) given 
above, fit) becomes linearly dependent upon b(t) and c(t), and we can thus, gauge it 
away -utilizing the remaining freedom in arbitrary constants, contained in Ap(t). 

Note: From now on, we drop the t-symbol -for time dependence-, from the various 
quantities; e.g., a stands for a(t). 

Inserting the form of 7 a/ 3 in equations (3.1b), we find that they vanish identically. 
We next consider, the following quantity q, which is scalar under a general linear mixing 
a a -> a a = S^a a , with S% G GL(3, 3?), 



q = c; u c^ KX r T i ua = - = 



2 7 2(bc-P) 



where 7, as usual, denotes the determinant of the matrix •y a/3 . Then, (2.7)., gives the 
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following non-zero components for the Ricci tensor R a p, and the Ricci scalar, R: 

Rn = -97n 

Rrs = qirs r, s = 2, 3 ( 3 - 3 ) 
R = q 

The (1,1) component of (3.1c), is an autonomous equation for the scale factor a: 



(-)'+a 2 = (3.4) 

(X 

with a first integral: 

(-) 2 + a 2 = uj = constant > (3.5) 

(X 

Using (3.2), (3.3) and (3.4), we get the equation for q: 

+ 3a 2 = (3.6) 



'9v , o 9 



q 

To obtain first integrals for (3.1c), let us define the new variables: 



7n = q 1/3 7n irs = q 1/3 irs 



(3.7) 



Then: 



Tyll = g l/3 7 ll = g-V3 7rs 



a 2 



7 = <fet(7«* ) = ^ 1/3 7 = jq~ 2/3 (3.8) 

It is straightforward to see that, with the use of (3.7), and (3.3), (3.4), (3.6), the spatial 
Einstein's equations (3.1c), translate into the following simple, integrable, Kasner-like, 
equations, -in terms of 7^: 

{r p i PP y = (3.9) 
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with first integrals: 



(3.10) 



where: 



( i \ 

e\ o o 



i9 Q - 
Vp - 



9 q 



\ 



a 7r 

Taking the trace of (3.10), we obtain -by means of (3 



while the (1,1) component of (3.10), gives: 



a 3 q 



(3.11) 



(3.12) 



The last two, imply that 6\ = $ s s = 9 + n, so finally, the matrix i? becomes: 

/ \ 

^ + 7T 



i9 a — 



9 q 



\ 







a it 



(3.13) 



Using the relation 7 = a 2 / (2g) -earlier mentioned- as well as (3.5) (3.6), (3.7) and (3.10), 
it is straightforward to see that the quadratic constraint equation (3.1a), becomes a 
relation among constants; that is: 



(3.14) 
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Integrating (3.5), we get the scale factor a: 



a(t) = a 



cosh(±A/cJt) 

^From this relation and (3.12), we conclude that: 

q-^ 3 a = a e^*, a > 



(3.15) 



(3.16) 



Now utilizing, in matrix notation, the relation: 7$ = -$ T 7 -which is the consistency 
requirement for (3.10)- and (3.16), we deduce that classical solutions exist, only for 
matrices with real eigenvalues -and thus diagonalizable. Since (2.16), (2.17) admit 
the solutions A^=constant, P a = 0, we can invoke a constant mixing of with a matrix 
of the form: 



A = 



V 











A3 


A|, 





/ 



and reduce it, to a diagonal form. Then, we are essentially led to the diagonal Taub's 
solution: 



7 22 = q^b = e et 



733 = <? 1/3 c = e nt 



(3.17) 



At this point, it is interesting to observe that, if for some reason, we had not invoked 
this diagonalizing A, and instead proceeded with the general we would had arrived, 
at a reducible form of the solutions with a non- vanishing 7 23 . However, this off-diagonal 
element, can be made to vanish through the action of the, previously mentioned, A. 
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Thus, we have shown that "gauge" transformations (2.14) -with (2.16) and (2.17), 
holding-, suffice to reduce the most general line element, for the Type II Bianchi Model, 
to the known Taub's metric. According to the theorem stated at the end of section 2, 
these transformations are, essentially, unique. We are now going to explicitly verify it 
-for the case at hand. 

A convenient way to proceed, is to start from Taub's form of the solution and ask 
ourselves, what is the form of the most general time-dependent automorphism A.p(t), 
which retains the form invariance of Einstein's equations (2.5) -written in the invariant 
basis. Since we know that A^ and A3, can be time-dependent, we focus on a time- 
dependent matrix A, of the form: 

/ \ 

£ 

A^= Oft Q 2 (3-18) 

where: g = Q1Q4 — Q2Q3, and all £>'s, are time-dependent. 

Consider the transformation, induced by this A|g, on r y^ ub -in matrix notation: 

7 = A T 7 M A (3.19) 

The linear constraint equations (3.1b), still imply N a = 0. As far as the time gauge 
fixing condition N = is concerned, we have: a/7 = \A\^/'y Taub , |A| > 0, and thus: 

N Taub dt Taub = ftdt^ 

dt\A\^ Taub = v / 1 Taub dt Tauh =>■ 
dtg 2 = dt Tauh 
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Since we wish for the transformation, to be a symmetry of (2.5), and we have secured 
that N a = 0, and selected N = a/7, the equation satisfied by j a p, would be exactly (3.1) 
and (3.2). Only the dot -defining the time derivative, with respect to Taub's time-, must 
be replaced by a prime: 

' = i= Q 2 (t Taub ) -j-^—r = Q 2 {t Taub ) X- (3.20) 

dt dt Taub e y ' y ' 

Defining the corresponding scale quantities 7 a/3 , -according to (3.7) and (3.8)- we must 
have the analogues of (3.10): 

Ti = n ( 3 -2i) 



Equation (3.2), reads: 



(Zr) ' ~ ^R = (3.22) 

7 

It also holds: 

( — )' - ^TaubRTaub = (3.23) 
iTaub 

Translating (3.22) in the t Taub - variable, and subtracting (3.23), we get: 

2(f ? 2 y + (£ 2 y— = o 

iTaub 

which, with the help of -f Taub = a 2 Paub /2q, (3.12) and (9\) Taub = ($ s s ) Taub , becomes: 

2(^)" + (^)-(2(^) T "*-^) = (3.24) 

Q 

The (1,1) component of (3.21) is: 

7n 
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where: 



9 «-l/3C? _-l/3-2. Taub 



and thus, that component reads: 



te 2 )' + e i W Tavb = & 



(3.25) 



Inserting the derivative of (3.25) into (3.24), we have: 



<✓)■£ = o 

which in conjunction with (3.6), implies £ 2 =constant. Without loss of generality, we can 
take g 2 =l. Henceforth, the time variable t, may -and will- denote Taub's time. 

It is thus left for us, to investigate the unimodular matrices: 

/ \ 

1 



Qi Q 2 



y £ 3 g A j 

with: 1 = QiQ A — Q2Q3, and all p's, are time-dependent. 

It can be proved that a convenient parametrization for this task, is: 



( 1 X 
1 








K 



where: 



a; = R r m LT 



33 



/ \ 

v(t) x(t) 



V 



o imo 

and R r m , are the symmetries of the Taub's metric, i.e. 
R T ^ Taub R = 7 Taub -in matrix notation-: 



( 1 ^ 

1 



R 



T3T 



V 



o 



/ 



pr 



R r m being: 

/ 

cos(£(t)) sin(^(t))e-( K -^/ 2 
sin(^(t))e( K -^*/ 2 cos(£(f)) 
where g(t), is an unspecified function of time, and k, /x, the eigenvalues of $ Taub . 
The system (3.21), gives the following equations for x(t) and ip(t): 



( 
\ 



\ 

/ 



(3.26a) 



if if if 



(3.26b) 



(3.26c) 



-2^ + /i = 7T - (7- 



Equation (3.25), for the choice f? 2 =l, gives ($ s s ^ Taub = and hence: 



(3.26d) 



it + 6 = k + jJ, 



(3.27) 
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It also implies, j u = 7n aub , or a(t) = a Taub (t), as well as, u = u Taub , or -through (3.14)-: 



2(9 + tt)' 2 + 7i9 — ga = 2(k + fi) 2 + k/j, k/j, = 9-n — go 



(3.28) 



Out of the 4 differential equations (3.26), only the first three, are independent -in view 
of (3.27). The solution to this system, for a ^ 0, is given by: 



— = ki 



A 3 c 4 e -Ai 



a(l + A 2 c 4 e- A *) 



A = K — fJ, 



(3.29) 



-from the Riccati (3.26b), where k\ = (n — 9 + X)/2a, is the constant special solution 
and: 



y? 2 = 



a 



(l + \ 2 c 4 e~ xt ) a>0 



A 2 c 2 



(3.30) 



Thus, it is easily seen that, (3.29) and (3.30) make the matrix L™, to be written in the 
form L™ = E™L™, where: 

cos(g(t)) sm(g(t))e' xt/2 ' 



V 



sin(g(t))e 



Xt/2 _ 



cos(g(t)) 



^ e 2 c^ (k ± - ^)e 2 Cy/a 



with (ei) 2 = (e 2 ) 2 — 1, (a, c) > and: 



tan(,(t)) = £H^ e -^ 



There are the special cases a = 0, or A = 0, which are easily seen, to fall into the 



previous case. 
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Thus, in all cases, there always exist matrices £ and L, such that the transformation 
matrix A^j, can be written as: 



\ 



V 



o 









pr v m T n 



J 



This concludes the verification of the Theorem stated at the end of the section 2, 
since indeed, R and S, have trivial action, on 7^ ufc . It is therefore, evident that the 
most general 7 Q( g, N(t) and N a (t), satisfying equations (2.5), are -in matrix notation: 



Imost general® = A T (t)^ (h(t))A(t) 



0104 -02 03 x{t) y(t) 



A 







V 



01 02 



03 04 



where, the £>'s are constant, and: 



'Jmost general \h(t) 



N a (t) = S a Jt)P?(h(t))h(t) 



pp(h(t)) = {p{t) 6lV ® - CM® - i 

' (0104 - 0203)h(t) ' (QiQ A - 0203)K t ) 



S = A-\t) 
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7 Tau6 (Mt)) 










e ( K +2 M )h(t) 



cosh(±y/tvh(t)) 



CO = 2(k + /i) 2 + Kjl 



where the fourth arbitrary function h(t), accounts for the time reparametrization covari- 
ance, i.e. permits us to depart from the time gauge fixing N = y/j. 



3.2 Bianchi Type V 

As it can be seen, from the results of section 2, concerning Type V, we can consider -with 
the usage of time-dependent A.I.D.'s-, the time-dependent part 7 Q| g(£), of the 3-metric, 
to be of the form: 

' a{t) b{t) 
lap{t) = b{t) c(t) 

v /(f) J 

with a(t)c(t) — b 2 (t) = f 2 (t). Again, as it happens for Type II, the form of the allowed 
transformation A^(t) is such that, one can not set b(t) = 0, a priori. Yet, we know -see 
(3.39) and (3.40) below- that the diagonal Joseph's metric, is the irreducible form of the 
most general Type V, solution. This puzzle, finds its resolution only on mass shell; after 
we have completely solved (3.1) with ^ a p{t) given above, b(t) becomes linearly dependent 
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upon a(t) and c(t), and we can thus, put it to zero -utilizing the remaining freedom in 
arbitrary constants, contained in A^j(i). 

Note: From now on, we drop the t-symbol -for time dependence-, from the various 
quantities; e.g., a stands for a(t). 

Inserting the form of 7 a/ 3 in equations (3.1b), we find that they vanish identically. 
We next define, the scalar -under a general linear mixing a a — > a a = S^a a , with 
S% E GL(3, 3?)- quantity q: 

The condition ac — b 2 = f 2 , now reads as: ac — b 2 = l/q 2 , or 

7 = ^ (3.31) 

Then, (2.7), gives: 

R aj3 = 2q>y afS 

(3.32) 

R = Qq 

The (3,3) component of (3.1c), gives an autonomous equation for the scalar quantity q: 

(f )■ + ^ = ( 3 - 33 ) 

with a first integral: 

(-) 2 = = constant (3.34) 

q q 2 
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Defining the scaled quantities: 



7a/3 = 97a/3 



ItI = i 



(3.35) 



and using (3.32), (3.33), equations (3.1c), are translated into the following form: 

(r p l P pY = (3.36) 



with first integrals: 



where: 



( 



(3.37) 



V 



e 


Q 


a 


-e 









\ 



The form of the matrix has been derived, using the form of 7 a(g and the property 
that |7| = 1. Using (3.31), (3.34) and (3.37), the quadratic constraint (3.1a), becomes a 
relation, among constants -as it was expected-, namely: 



3co = 6 2 + QO 



(3.38) 



The property |t~| = 1, together with the consistency requirement -in matrix notation- 
y# = ■d T7 y, which follows from (3.37), enables us to conclude that classical solutions, 
exist only for those values of the parameters, 9, g, a, for which is diagonalizable, i.e. 
when 9 2 + ga > 0. 
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Since the matrices A^, of the form: 



A /3 



/ \ 

Qi Q2 
Q3 Qa 



v o 1, 

along with P a = 0, constitute the remaining gauge freedom, we can invoke such A% to 
diagonalize fig, and -at the same time- retain the shift, zero -see (2.14). Now with a 
diagonal equations (3.37), essentially imply that 7 Q/3 , is diagonal too. 
A further integration of (3.34), yields: 



/(*) 



^sinh(± v /cJt) u > 



(3.39) 



±2t 



and thus, we are laid to the well known Joseph's solution -through complete integration 
of (3.37), for the diagonal case: 



7n = qa = e 



7 2 2 = qc = e 
3u = X 2 > 



(3.40) 



or the Milnor's solution ||18| ) when uo = -with the corresponding q. 

Once again, it is interesting to observe that if, for some reason, we do not invoke 
this diagonalizing A? and, instead, proceed with the general we arrive at a reducible 
form of the solution, which contains a non- vanishing j 12 . However, this off-diagonal ele- 
ment, can be made to vanish through the action of the -previously mentioned- constant 
automorphism. 
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Thus, we have shown, that the "gauge" transformations (2.14), -with (2.16) and 
(2.17), holding- suffice to reduce the most general line element for the Type V Bianchi 
Model, to the known Joseph's metric, as predicted from the theorem, stated at the end 
of section 2. As we have done for the Type II case, we are now going to explicitly 
verify that these transformations, are essentially unique. To this end, let us consider 
the most general time- dependent automorphism, complementary to the time-dependent 
automorphism, described in section 2 -for the Type V, case. 



A /3 



/ A(t) B(t) ^ 
C(t) F(t) 



v 1, 



(3.41) 



with A(t)F(t) — B(t)C{t) = 1. The action of such automorphism on 7^ seph , is -in matrix 
notation: 

j = A T 7 Josep/l A 

If we insert 7 aj a, in the linear constraint equations (3.1b), we learn that N a , are also 
zero and, since, \jap\ — IM^lap^l = ITo/T^'I, we conclude that we are in the same 
temporal, as well as spatial, gauge. Therefore, 7 Q| g, will also satisfy equations (3.1c). 
Since Ajg, is an automorphism, it is a symmetry of q and thus, if we define the scaled 
quantities: 

lap = Wad 

they must satisfy, the relations analogous to (3.37): 

Tl pP = n (3-42) 
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where: 



,9" 



9 g 
a -6 



y j 



while, 7„^ ep/l , satisfies the relations: 



( 7 T P ) Joseph {ipp) Joseph — fs) Joseph 



where: 



A 



in) 



/3 ) Joseph 



-A 



V 



/ 





By virtue of (3.34), -and since q, is invariant-, we get that uo = UJ Jose P h ) i. e . 



9 2 + ga = A 2 



(3.43) 



In order to proceed with the integration of (3.42), it is convenient, to parametrize 
in (3.41), as follows: 



/ 



A /3 



\ 



M 



1 



with A' = RLL™, where R r m is: 



V 



e- A '/ 2 ^ ( cos(g(t)) sm(g(t)) 
e xt ' 2 J y -sm(g(t)) cos(g(t)) 
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e xt/2 o 
e- xt ' 2 



I 



i.e. the symmetries of the Joseph's metric; -in matrix notation 

R T 1 Joseph R = ^Joseph &nd jja ig . 

/ \ 

<p(t) r(t) 
K l/ip(t) } 

The system (3.42), gives the following differential equations for <p(t) and r(t): 



2^ + A = 9 + a- 



(3.44a) 



Tip — rip = gip 2 — 29ipr — or 1 



(3.44b) 



e 2Xt (ip-rp) = — 



ip* 



(3.44c) 



The solution to this system, for a ^ 0, leads to incompatibility of the form ip 2 = — e 2 , e 
a function of time. 
For a = 0, we get: 

ip(t) = cie^* 

(3.45) 

T(t)= Cl f e e^ 

with ci > 0, and -from (3.43), for the case at hand-, 9 = ±A. The case 9 = A, trivially 



gives, a constant matrix 



/ 



Ci Ci£ 

1/ci 



while, the case 6 1 = —A, gives 



( 
\ 



e~ A * 
e xt 



l/ci 

Ci -Ci^ 
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Since the first matrix in the product, is a symmetry of (7^) Joseph, we again conclude 
that, the non-trivial action of A^j, on (7^), Joseph, is tantamount to the action of a constant 
matrix in accordance to the theorem of section 2. 

Finally, the most general line element (7^/3, N, N a ) satisfying Einstein's equations 
(2.5), is thus given, -in matrix notation by: 

l m ost general® = A^h^tMO) A(t) 



Ql P(t) g 2 P(t) x(t) 



A 



QsP(t) QiP(t) y(t) 
1 

where the p's are constant, subject to the condition — Q2Q3 = 1 and: 



V 



/ 



N(t) — J \^ m ost general\h(t) 



N a {t) = S%pP{h{t)) 



Pf(h(t)) = {x(t)(\n j^YMmn |^)-, (In -^) } 



r Joseph (h(t)) 



S = A- 1 

( e Xh(t) 

q 










\h(t) 





q 







./ 



\ 
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1 



q(h{t)) 



(3.46) 



f(h(t)) 



±2h(t) 



u = 



3u = A 2 



where the fourth arbitrary function h(t), accounts for the time 



reparametrization covariance. 



4 Discussion 



In this work, we present an approach to the problem of solving Einstein's equations, for 
the case of a generic Bianchi-Type spatially homogeneous spacetime. The approach is not 
plagued by the fragmentation characterizing the major part of the existing rich literature 
-which is inherited by the diversity of the various simplifying ansatzen, employed in 
each case. The key notion for avoiding this fragmentation, is that of a Time-Dependent 
Automorphism Inducing Diffeomorphism; that is, a general coordinate transformation 
(2.9), mixing space and time coordinates, whose action on the line-element of a Bianchi 
Geometry, is described by relations (2.14) -viewed as "gauge" transformation laws for 
the dependent variables 7 Q| g(i), N(t) and N a (t). The investigation for the existence 
of such G.C.T.'s, leads to the necessary and sufficient conditions (2.16), (2.17); hence 
the name Time-Dependent A.I.D.'s. In each and every Bianchi Type, these conditions 
possess a non-empty set of solutions containing precisely three arbitrary functions of 
time. A choice of these arbitrary functions, amounts exactly to a choice for the three 
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spatial coordinates. Thus, the possibility is offered for simplifying Einstein equations, 
-through a simplification of j a @, N a , N-, without running the risk of loss of generality 
or any sort of incompatibility. 

Of course, the possible simplifications differ from one Bianchi Type to another; even 
within the same Bianchi Type, there are many possible simplifications -since one, can 
use the three arbitrary functions at will. This kinematical freedom, when combined to 
the dynamical information -furnished by the linear constraint equations-, considerably 
simplifies the form of the line-element and thus of Einstein's equations, as well. A useful, 
in our opinion, irreducible form of the line-element for each Bianchi Type, is given at 
the balance of section 2. 

A statement that applies to all Types is that, using two of the three arbitrary functions, 
the scale-factor matrix "fapif) can always -a priori; i.e. before solving any classical 



equations of motion- be put into a so called "symmetric" [16| form, i.e. 713 = 723 = 0. 
This applies also for Type II, if we take instead of the standard form for the structure 
constants (C23 = —C\ 2 = 1, all other vanish) the equivalent version Cf 2 = — Cfi = 1; an 
other vanish. If this " symmetric" form, is then substituted into the linear equations, and 
the third arbitrariness is used, considerable restrictions among N°"s and the remaining 
7 a/ g's are obtained, as presented in detail at the end of section 2. Furthermore, with the 
help of the essential arbitrary constants in Ajg, we can diagonalize 7 aj g(t), on mass-shell. 
For all Bianchi Types, the shift vector N a , can always be set to zero -with the help 
of Time-Dependent A.I.D.'s, and the linear equations. One could of course, rely on 
the well-known existence of Gauss-normal coordinates || , and argue that this should be 
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true. However, in this work, the explicit realization of this fact is presented; what is more 
important, is that the vanishing of N a , is accomplished without spoiling manifest spatial 
homogeneity. The interplay between line-elements with and without shift, established 
through Time-Dependent A.I.D.'s -see (2.14b)-, raises the need to reexamine the set of 
existing solutions -with respect to physical equivalence, among each other. In particular, 



many tilded and untilded fluid solutions fL8| , may proven to be G.C.T. related -and thus 
physically indistinguishable. 

Except of the three arbitrary functions of time, of considerable importance, are also 
the (non absorbable in the shift) arbitrary constants, appearing in the solutions to (2.16) 
and (2.17). The number of these constants, varies for different Bianchi Types. The very 
interesting fact, is that when this number is subtracted from the number of constants, 
given by Peano's theorem, -after the freedom in arbitrary functions of time, has been 
fully exhausted-, the resulting number of the -finally- remaining constants, equals, for 
each and every Bianchi Type, to the number of expected essential constants -see [ |19j] , 



p. 211. This, permits us to conclude that the gauge symmetry transformations (2.14) 
-with (2.16) and (2.17), holding- are, essentially, unique. It is also noteworthing, that 
the existence of these constant parameters, helps to rectify a defect from which, the 
previous approach of Jantzen, is suffering; that of an uneven passage, from the lower to 
the higher Bianchi Types, owing to the change of the dimension of the invoked symmetry 
group [15]; indeed, the arbitrary functions of time are thus varying with dim[SAut(G)], 



from 8 (Type I), to 5 (Type II and V), to 3 (higher Types). This situation, is rather 
unsatisfactory, since we know that the independent or dynamical degrees of freedom 
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for the gravitational field, are 2 -per space point. Thus, in cosmology, we expect 2 
independent functions of time -irrespective of Bianchi Type. 

In contrast to this state of affairs, the solutions to (2.16) and (2.17), contain exactly 3 
arbitrary functions of time, which together with the arbitrary function -owing to the time 
reparametrization covariance of equations (2.5)-, leave us with 6(7 a/ g) — 4 = 2 arbitrary 
functions, in all Bianchi Types. The required sensitivity, of the method, to the particular 
isometry group, is represented by the extra constant parameters -as explained. 

It is in this remarkable way, that General Relativity manages to encode the memory of 
spatial G.C.T. covariance, in the set of the reduced equations (2.5), where only functions 
of time and their derivatives appear. This encoding also persists in the actions -when 
these actions exist-, and leads to important grouping of 7 a/ g's, into the three solutions: 
x 1 = C'^C , p (J 7 Atp 7 y<T 7a/3, x 2 = C^ & Cl a ^ v , x 3 = 7 of the quantum linear constraints 



fl3| , When a truly scalar Hamiltonian exists [13], the wavefunction depends only 
on the g*'s: 

q Vl 2 7 4 ^ 7 7 7a/3 ' q ^7 



which completely and irreducibly, determine a spatial three-geometry. 

To summarize, the system (2.5), admits solutions containing in each and every 
Bianchi Type, exactly four unspecified functions of time. One, corresponds to the free- 
dom of changing the time coordinate; three, correspond to the freedom of changing the 
spatial coordinates via Time-Dependent A.I.D.'s. The action of such a transformation on 
the line-element, and on the system of equations (2.5), is described by relations (2.14), 
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(2.18). Thus, one does not actually need to calculate the simplifying G.C.T.'s; one sim- 
ply uses (2.14), simplifies the equations, solves them completely, and finally inverts the 
transformation thereby obtaining the entire space of solutions. It is in this sense, that 
the closed form of the line elements presented in section 3, exhaust the space of classical 
solutions -for the case of Bianchi types II and V. 
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The number of arbitrary constants appearing in general solution for each Bianchi 
Type -vacuum model-, is given in the following table -depicted in the first of [H\] , pp. 
211: 



Bianchi Type 


# of the essential constants 


I 


1 


II 


2 


VI Q , vii 


3 


VIII, IX 


4 


IV 


3 


V 


1 


VI h (h -1/9) 


3 


W-1/9 


4 


vn h 


3 
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A Appendix 



In |[22| , the following Theorem, is given: 

"Let two symmetric forms A and B, be given, on a n- dimensional linear vector space V. 
If one of them -say A- is non singular, then there is a base in V in which both A and 
B, are diagonal, if and only if, the mapping A~ 1 B, possesses n-real eigenvalues. " 

Thus, if we take the pair j a p, r] a p, it suffices to prove that T] a9r y e /3, has n-real eigenval- 
ues. In what follows, for the sake of completeness, we give a proof of the entire statement 
that a positive definite matrix 7^, can be diagonalized via the Lorentz group. Theorem 
Let 7 be a positive definite n x n real matrix. Then, there exists a Lorentz matrix A, 
such that: 

A T 7 A = A (A.l) 

where A a diagonal matrix. 

Note: Since A T = r]A~ 1 i], where rj is the Minkowski metric, (A.l) may be written as 

A~VyA = rjA (A.2) 

In order to prove (A.2) it is useful to write it equivalently using the notation employed 
with linear mappings. To do that, we consider an n-dimensional real linear space V with 
basis (ei, e 2 , . . . , e n ). The scalar product in this space is defined as < , > : V xV ^ 
with < e a , ep >= rj a p. The matrix 777 defines a mapping / : V — > V through the 
relation: 

n 

/Oa) = ^{vrdapep 

(3=1 
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The following will prove useful later on: 
1) If M C V then V = M © M 1 - [|| . 



2) A mapping / : V — > V is called self-dual, if 

< /(x), y >=< x, f(y) > for every x, y E V. We may prove that the mapping / 
defined through the matrix 777 is self-dual. Indeed: 

< f(x), y >=< y, f(x) >= y T r]T]-fx = y T jx 
< x, f(y) >= x T r)T)jy = x T 'jy = y T ^x 

=> < f(x), y >=< x, f(x) > 

3) If M C V is an invariant subspace of V with respect to a self-dual mapping / then 
M x is also an invariant subspace of V. Indeed, let b E M x and m E M. Since M is an 
invariant subspace, it follows that: 

f(m)EM^> <f(m),b>=0^ < m, f(b) >= V m E M 

=> f(b) E M x 

Equation (A. 2) states the fact that there exists an orthonormal basis of V consisting of 
the eigenvalues of /. If (A. 2) holds then the non-vanishing elements of the real diagonal 
matrix 77A will be eigenvalues of 777. Thus, we have to prove that the eigenvalues of 777 
are real. Indeed, the following theorem holds: 
Theorem 

If 7 is a positive definite symmetric matrix, then 777 has real eigenvalues. 
Proof 
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Let A = a + (3j, /3 ^ a complex eigenvalue of 777 and u 7^ the corresponding complex 

right eigenvector. Since 77 is invertible, there exists a v = x + yj, x, y, G 9ft n such that 
u — rjv. We have: 

777W = \u 77777-1; = A777J 

77777a; = arjx — j3r]y (A. 3) 

VIVV = ar W + (A-4) 
Equations (A. 3), (A. 4) imply respectively: 

y T r]-fr]x = a<y, x> -(3 <y, y> 

x T r\ r )r\y = a < x, y > +(3 < x, x > 
The last two equations have their left-hand sides equal (since 77777 is symmetric), hence: 

(3(< x, x > + < y, y >) = =^< y, y >= - < x, x > (A. 5) 
Since 7 is positive definite, 77777 is positive definite as well. Then: 

x T r]^jr]x > a < x, x > —(3 < x, y > > (A. 6) 



7/ T 77777|/ > ^=^a < y, y > +(3 < y, x > > ^=>^ 



a < x, x > —(3 < x, y > < 
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(A.7) 



^From (A.6), (A. 7) we get: 

a < x, x >= (3 < x, y > (A. 8) 

Through (A. 5), (A.8), equations (A. 3), (A.4) imply: 

x T r}^r\x = 

y T vivv = o 

and, since 77777 is positive definite we conclude that x = and y — 0, i.e. u — 0, 
contradicting our initial assumption u 7^ 0. Therefore /3 has to vanish and thus we have 
proved the reality of A. 

For the eigenvectors of 777, we can prove that they have a non-zero norm. Indeed, let x 
be an eigenvector of 777, i.e. 

777a; = \x 7X = \l]X =>- x T ^x = \x T 7]x = A < x, x > 

Since 7 is positive definite and x ^ we have x T 7x > 0, so that < x, x >^ 0. 

We are now in position to prove a spectral theorem for a mapping / with real eigenvalues. 

Theorem 

If / : V — > V is a self-dual mapping with real eigenvalues, then V has an orthonormal 

basis consisting of the eigenvectors of /. 

Proof 

Let A be an eigenvalue of /, u the corresponding eigenvector and M = [u] the one- 
dimensional subspace spanned by u. Obviously, M is an invariant subspace of V with 
respect to /. 
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According to 1), we have V = M © M L . As implied by 2) and 3), M L is also an 
invariant subspace and thus / induces a self-dual mapping onto M L . Hence, we can 
apply induction and show that: 

V = Mi © M 2 © . . . © M n 

where the M a are one-dimensional invariant subspaces orthogonal to each other. Since 
u a is an eigenvector of 777, it holds that < u, u >^ 0, as proved above. We can thus 
promote the orthogonal basis to an orthonormal set (ui, u 2 , ■ ■ ■ , u n ). The transformation 
connecting this orthonormal basis to the initial orthonormal basis (e 1 , e 2 , . . . , e n ) is the 
matrix A sought for in the first theorem, relations (A. 2) and (A.l). 
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